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Abstract. In this paper, the author propose a numerical method to compute the solution of the 
Cauchy problem: wt — {w™'Wx)x — w^, the initial condition is a nonnegative function with compact 
support, m>0, 1 < p < m + 1. The problem is split in two parts: A hyperbolic term solved by 
using the Hopf and Lax formula and a parabolic term solved by a backward linearized Euler method 
in time and a finite element method in space. Estimates of the numerical solution are obtained and 
it is proved that any numerical solution is unbounded. 



1. Introduction 

In this paper, we study a numerical method to compute the solution of the Cauchy problem : 

, . Wt — {w'^Wx)x = w^, t > 0, X G M 

^ ^ w{x, 0) = Wo{x) > 0, X e M 

Wo is a function with compact support, m > , 1 < p < m + 1. 

Samarskii et al jlO|, see also [I], [2], [2], [6], [5], [1],|9] have obtained theoretical results on this 
problem. In the case 1 < p < m + 1, the numerical solution blows up in finite time and there is no 
localization (HS-regime) that is u{t,x) — ^oo in R if t — >Tq. 

A numerical method to solve fll.ip has been proposed in the case of S-regime (p = m + 1) in [7] 
and in the case of LS-regime (p > m + 1) [Hj. If we denote ujl = {x G ]R/m(Tq", x) = oo] , in the 
case of S-regime, L* = mes{ujL) is positive while in the case of LS-regime , L* = 0. The problem 
is solved by using a splitting method; for that, it is more convenient to work with the function 
u = w^. Problem (11.11) may be written: 



^ ■ ^ u{x,0) = Uo{x) = w^{x), x G M 



with q = m > 0, q < 1. 

This problem is split in two parts: a hyperbolic problem which will be solved exactly at the nodes 
at each time step and allows the extension of the domain and a parabolic problem which will be 
solved by a backward linearized Euler method which allows the blow up of the solution. 
In [7j, [S], the convergence of the scheme has been proved in the cases q = 1 and g > 1. It has also 
been proved that for q = 1, the numerical solution blows up in finite time for any initial condition 
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and that its support remains bounded if the initial condition is smaller than a self-similar solution. 
In the case, 1 < q < any numerical solution is unbounded while for q > {p > m + 3) if 

the initial condition is sufficiently small, a global solution exists and if q > ^^^^ for large initial 
condition, the solution blows up in a finite time. We observe numerically that in any case, the 
unbounded solution is strictly localized and blows up in one point and that for q = ^^^^ also, if the 
initial condition is sufficiently small, a global solution exists. 
Here, we generalize this method to the case q < 1. 

An outline of the paper is as follows: 

In Section 2, we present the numerical scheme. In Section 3, we obtain estimates of the ap- 
proximate solution and of its derivative in space. In Section 4, we prove that if g < 1, (p < m + 1)) 
any numerical solution is unbounded 

2. Definition of the numerical solution 
In order to solve problem fll.2p . we separate it in two parts: a hyperbolic problem 

(2.1) ut - —(u^f = 0, X e M, t>0 

m 

and a parabolic problem 

(2.2) Ut - uu^x = mu'^^^, X eR, t>0 



We denote by Atn the time increment between the time levels t„ and tn+i, n > and by 
the approximate solution at the time level t„. This solution will be in a finite-dimensional space 
which will be defined below. 

Without loss of generality, we can assume that the initial condition is a continous function with 
a symmetric compact support [—sq, sq]. Let iV G N, the space step h is defined by /i = ^; we note 
Xi = ih, i E Z, Jj = {xi-i, Xi) and we define the finite dimensional space by 



(2.3) 



= {(PHe C°(M)/0^(x) = 0, X ^] - So, +Sol (j)h\i, e Pi, 2 = -iV + 1, N} 



For (ph € Vl^, we note (pi = (j)h{xi), i E Z and for any function v E C°(M) with compact support 
[—So, +so], we define its interpolate by vr^w E V^f and vr^f (xj) = v{xi), i E Z. 

The support of the solution will be denoted [— s", s"] and will be computed at each time 
level by solving (12. ip . 
We denote 



'si' 












_h _ 


_h _ 



(N'l - l)h, h'^ 



{Nl - l)h 
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(for x G M, [x] denotes the greatest integer less than x), so we get: h < h^^h^ < 2h. 

We then define the finite-dimensional space V/^^ by: 

e C^{R)/Mx) =0, x^]- si, sil 
V^={ <Ph\uePi,i = -N^ + 2,Nl-l 

and denote by vr^f the Lagrange interpolate in of a function v G C°(]R) with a compact support 
in [-si, s\]. 

The solution m^"*"^ at the time level tn+i is computed in two steps: knowing u^, we compute first 

an approximate solution of (12. ip which will be denoted ^ ; then starting with this intermediate 
value, we compute an approximate solution of (12.21) at time level tn+i- 
Then if S is the semi-group operator associated with (12.11) . we define 

= 5(At„)<. 

The support of this function will be denoted [— s1^^, s"^^] and the interpolate of m^+s in V^^'^ will 

be denoted u!^^ = 'k'^'^^u^^'^ . 

Then starting with u^'^ , the function is obtained by solving (12.21) with a backward lin- 
earized Euler method in time and a Pi-finite element method with numerical integration in space. 

This function has the same support as u^y^^ ■ 



2.1. Computation of the solution of the hyperbohc problem. 

The hyperbolic problem is independent of g; we use the same method as in [7] for the case q = l.It 
is not necessary in this case to use P2-interpolation on the last interval since there no localization. 
We use the Hopf and Lax formula which gives explicitely the solution to (12.11) with the starting 
data at the time level t = tn- Here, we simply recall the results obtained in [7] 

We define the piecewise constant functionwj^ by 

^ t »-l J. ^ _^n j^2< I <Nl- I 



U 



V 



n 



- — on (-si, X-N^+i] 



Let us denote r„ = \\v\\^ = \\v\\i^sm) , s > 0. 
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Proposition 2.1. // the following stability condition 



(2.4) rMoo<'^ 
is satisfied, then the solution u^^^ of l{2. 1\) is defined by 

uT'^ = < + ^ (max(0, t;r+i))' , -iV!5 + 1 < ^ < iV^ - 1 



At, 



s7 - 



m 



n+l _ „n I ^^I!:,," 

m - 



Ifm+^ = Nl + l, we get 



and we get analogous formula at the other end of the support. 



2.2. Computation of the parabolic problem. 

The approximate solution at tn+i is now obtained by solving problem (12.20 . 
We introduce the approximate scalar product on Vj^~^^ : Wcphyiph ^ ^/T^^? 

i=N"+^-2 

1 1 



We define m^"*" as the solution of the following problem: 




The second member of f l2.2p is splitted in two parts in order to obtain the L°°-estimate of 
This equation may be written: 
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Tl-\- — I 

^ + m(l - q)Atn [u 



1\ 9+1 



t+i 



.1 \ "? 



1 - mqAtn ( ) ) 



-u. 



h "^r^-i \hl+^'^K^'-^ h + h 



n+l 



1 — mqAtn ( M 



n+l 



n+l 



n+i X 



-Af!!+Hi 



We get immediately the result: 



fi-\-— 

Proposition 2.2. If the hypotheses o f proposition l27B are satisfied and if U/^ ^ satisfies 



(2.6) mqAt 
then the solution u^'^ of Ii2. ^j) is unique and nonnegative. 



< 1 



2.3. Properties of the scheme. 

Lemma 2.3. If the hypotheses of proposition l27B are satisfied, then the following estimate holds: 



(2.7) 



I ft 1 1 oo — 



\K\Li^ + mil-q)At^\ \u 
l-mqAtn\\ul\\l 



n\\q 
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Proof: We get immediately from the Hopf and Lax formula that 



io the index such that m"^^ 

^0 



n+l I 



we get from O, O that 



n+4 



Ui. 



< \\ul\\^. If we denote 



1 \ 9 



1 — mqAtn 



rill? 
oo 



which proves the lemma. 

We deduce the following theorem: 



Theorem 2.4. Under the hypotheses of proposition \2. 4 the numerical solution exists at least until 
the time 



(2.8) 



and the following estimate holds: 



1 



II nil? 



(2.9) 



l<lloo< 



Proof: This result is proved recurently. It is true for n = 0. If we suppose that we have estimate (12.91) 
at the time level t„ , we get from (12. 7p , at the time 

tn+l: 



,n+l II 
''h Woo 



l+m(l-g)At„Kr^ 



1 — mqAtn 



n\\1 
oo 



or 



I "■+111 ^11 II 

kit < \\Uf,\\ 

I fi' 1 1 oo — II 1 1 oo 



l-mqt^\\uUl + m{l-q)At„,\\uUl 
(1 - mqtn+i IKIIL -) (1 - mqtn iKIlL) 
The inequality (12.91) will be satisfied at the time if : 



(l — mqt 



n+l \\t^h 



0\\1 



+ mAt^\\ul\\l) 



< (l - mqtn \\Uf^ 



'^^')Hl-mqt.^,\\ul\\ir'^ 



I OO'' 



By using the Taylor formula, we get: 

(l - mqtn \\ul\\U' - {l-mqtn+i > mAtn \\ul\\l^ (l - mgt„+i ^ 
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and the inequality f l2.3p is satisfied: 



Lemma 2.5. Under the hypothesis of Proposition({27B), we have the estimate: 

\K\L < c 

for tn < T < Ti where C is a constant depending on T and uq. 



Proof: We have the inequahty( [7] ): 
It remains to estimate fr"*"^ 

II n, 1 1 oo 

From fl2.2p . we deduce the following equation satisfied by f^^^: 



i-mqAt^[ur'^y)vr'+^ 



,n+l / „,"+ 



At,- 



mq- 



h 



n+i 



l\1 



i\ 9 



-(C7) )<„Y + ^(i-?) 



Atr 



n+4 



IN 9 + 1 



IN <J+1 



71+1 



'i<i<N 



n+l 



and we have analogous inequalities for i = N^^'^ - 1, Nl+\ -AT^+i + 2, -AT^+i + 1. 



By using fl2.2p for i — 1, we can replace u^^^ in the second member by its expression in function of 
the values of u^^'^ and f^^^ and we get: 



l-mqAtn(u''''~A')v:+' + 



Atn 



( 



1 n+ 1 ^ ~ mqAtn ( ( M,. 

v7 ' " — v7 1 u 



n+4 



1 \ q 



1 — mqAtn I u. 



n+4 



1 \ 9 



n+^ , Atn I I 



1\ 9 



1 1 + m(l - q)Atn [ul^^ 



i\ 9 

2 



1 - mqAtn 



-|-m(l — g) 



At, 
h 



n n 



IN 9+1 



- u 



IN g+1 



"i-i 
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Let io the index such that 



n+l 
«0 



l^^'''^lloo ■ -^^^^ previous equahty, we get: 



(l-mgAt„,Krj|</i| 



< 



n+4 
«0 



+ m(l - g)^— 



2 



«0 



- U 



n+4 



+mq 



U. 



„+i 1 +m(l - g)At„ IIm 



n+l II 
h Woo 



1 — mqAtn \\u 



n 

h II oo 



We easily obtain: 



and 



u. 



«0 



iO-1 



<(g+i) IK 



n||ij 

oo 



1X9 



Mi 



1\ <? 



^.0 



^ < 



«0 



Ml 



|9 

I oo 



and we get: 



(l-mgAt„K||^) \W'\\< 



a l + m(l-g)At„||<+^||' 
1 + m(l - g^) At„ l|<r^ + mgAt„ il^^At^Kf 

By (12.91) . we get easily there exist a positive constant C depending on m, g, T, such that 

Hence for tn < T < Ti, we get ^ C*. 

Lemma 2.6. Under the hypotheses of Proposition Ii2.^) . we have the estimate: 



(2.10) 



l<lli<C 



/or tn < T < Ti where C is a constant depending on T and uq. 



Proof: From the properties of the semigroup operator S [7J, we get: 



n+4 



< Wv^W-^, and by 



using the equations satisfied by f"'*'^, — A^""*'^ < i < NT^'^, we obtain: 
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1 — mqAtn 



n+4 



U^. 



,n+l|| 
111 



< 



and we immediately deduce the estimate (12.101) . 

In this case, since g < 1, the variation of is not bounded. 



3. Blow-up of the solution 



In this part, we prove that for g < 1 the solution blows up in finite time. 



3.1. Construction of unbounded solutions. 

Define the function 



d{x) 



1 — %, |x| < a 
0, Ixl > a 



rO 



We note 6h = fChd- If the initial condition is uj^ = -\Ofi{^ ) with ^' 

^ -OkiCn), with 



possible to choose A and a in such a manner that u'^{x) > 



q_t , we prove that it is 



1-9 



(T-t„)9 



= X {T — tn)) '^"^ and then the numerical solution blows up in finite time. 

q-1 



— Gh{in)- 



We denote = (T - t„) 2g 

{T-t„)i 

The support of u'^ is] — a(n, aCn[; its lenth is increasing with the time. 
If < u^, we get < 



n+i 



The support of i^'^ is s^^'] with s^^' < s 

Since is a symmetric function, it is sufficient to study the case x > 0,{i > 0). 
We get for i > 0, 



< S 



n+1 



At. 



m 



for z such that Xi < a(n since the function is decreasing for x > and we have: 



h 



{T - Q 



h 



with 9^ = d{Cl). 
Hence, we obtain: 
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ur = — I + 



with r , = J, = I (e + e-i)- 

Then ul+^ will be a subsolution of if 




< ^ + m(l - q)Atn iu 



By using the equality -.^ {26^+' - 6^+' - 6^^,') = 
this inequality reduces after simplifications to 



+m(l - q)Atn (^u 



Noting that i = + ^^'^ i = 0? + V? with bfl < we get: 



* V * V ma^T — tnJ ma^ T — tn 

and the inequality (13. ip becomes: 



^^^^ (i - (1 -„..) fC'l'T" 



At„ / „+i\9+i 
+m(l-g)^(^M. 2j 



By using the equality 9"''^^ = 6^^^ + 1 — 7!^, this inequality reduces to: 

gn,^ I 4A r-t„+i ^ 2A / 4A At. 



ma^ T — t„ V T — t 
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^ 4A r - t.+i / _ 2A At„ _ A _ /C+i _ ^ T-t^ 
" ma^ T -tn V a?T - t^+i ' / V Cn / 

If we denote Hn = — -o ~ since u""*"^ > — ^-^^j" (1 — l^'n), the preceding inequality will be 

ma T — tn {T-t„)i 

satisfied if: 



9 ' 9 r^" \ ' 9 / / — \ 9 f^" \ ' 9//\ '«/ \ A9 / Aj. 



/•2 / a2 ^2 



+mgA^^ (C)" (1 - C7i^ + 1 - .2 



bn+1 



(3.1) + m(l - g)A'^^^^ (^D^^^ (1 - -"n)^^' 

From the stability condition, we get: 

h h 

hence for h sufficiently small, we get: ^ — yU„ ( 1 H — > 0, and Irjl'l < 5 

ma \ a / 

So the inequality (13. ip will be satisfied if : 



+mgA'' (1 - ^n)-^ {O'lf + ^ - 1 + m(l - g)A^ (1 - /i„ 



Since ^" G (0, 1), we introduce the function $„(?/) defined on (0, 1) by : 



$„,(y) = mA^ (1 - y'^+i g + (1 - - 



T — f 



ma^ \^ cfi J J V Cn J At 
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4A 2A / 2A\ 



or $„(y) = AAV+' + C-% 
with 



A = m{l-^^r,Y(q+{l-q){l-^i^)^ 



T-t 



n 



4A 2A 2A 



5 = 1 + — - + — + — 



^ , 4A / 2A\\,^ ^, A T-t. 

^= zr-o-z^n 1 + (1-5) 



ma^ \ O? J J V Cn / ^'^n 

A sufficient condition to satisfy (13.11) is: $n(2/) > 0, ?/ G (0, 1). 
We have $„(0) = C, 
hence we get 



^ ma^ \ 0? J J V Cn ) l^t 

But since < g < 1, we get: 

< ^ _ 1< ^ " ^ 



Cn Q T — tn+l Cn+1 

and 

\ma'^ \ J J q 

So, if the quantity 4- is such that — ^ > ^, if is sufficiently small, we get C > 0. 

" ma g 

Let us define yo = ^, ?/o G (0, 1) and if y < yo, we obtain > 0; if < 2/ < 1, we obtain: 

$n(2/)> AA%^+^-i?(l-l/o) 

This quantity will be positive if : A^ > ^~q+\ ■ 
Further, we have: 



A>m{l-5Y[q + {l-q){l-5)[l- 



T-U 



The solution at the time level tn+i exists if mqAtn IIm^II^ < 1, that is — — — < — : 

and if A > Ao > — we get: A > mq{l — 5y 

(mq) 1 

and ^n{y) will be positive if: 
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(3.2) 



A" > 



mg(l — 5y \q 



1 2A 4A , 
- + ^ + — ^6 



ma 



1 + 



4A 2A 

ma a 



4A 



ma 



'l-6)-6 1 



2A 



ma 



The second member is a function of hence if \ is fixed such that 



4A 



> 



1-g 



ma q 

the inequahty fl3.2p will be satisfied if A is large enough and h sufficiently small. 
Hence if the initial condition satisfies uohix) > -\Oh{C,o), A and a satisfying (13.21) . the solution blows 

up in finite time. 



J'q 



3.2. Blow-up of the solution. 

Theorem 3.1. Let < q < l.The solution of problem blows up in finite time. 



Proof:Since Uo{x) ^ 0, there exists p > 0, e > 0, Xq such that uq{x) > e > for x such that 
|x — xo| < p. So we can choose T large enough such that ^ < e and f_q < p. 

We have Uq{x) > -\0h (^^-^j then the solution blows up at time Tq < T* where 
T* = Max ( (^)'' 



_2q_ 
1-q 



In Figl, we present the evolution of an initial condition uq for m = 1, p = 1.5. The solution blows 
up in finite time. 




Figure 1. blow-up for m=l, p=l,5 
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